function p_vec = prob_vector(data, sample_freq, first_freq, second_freq)

% function p_vec = prob_vector(data, sample_freq, first_freq, second_freq)

p_vec = zeros(length(data),1);

for k = 1:length(data)
    G = parameter_matrix(length(data), k, sample_freq, first_freq, second_freq);
    p = marginal_posterior_probability(data, G);
    p_vec(k) = p;
end